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Using the APM cluster data we investigate whether the dynamical status of clusters is related 
to the large-scale structure of the Universe. Due to difficulties in determining unambiguously 
the dynamical status of clusters in optical data, our substructure identification method has 
been calibrated using jointly ROSAT X-ray and optical data for a subsample of 22 Abell 
clusters. We find that cluster substructure is strongly correlated with the tendency of clusters 
to be aligned with their nearest neighbour and in general with the nearby clusters that belong 
to the same supercluster. Furthermore, dynamically young clusters are more clustered than the 
overall cluster population. These are strong indications that cluster develop in a hierarchical 
fashion by anisotropy merging along the large-scale filamentary superclusters within which 
they are embedded. 



1 Overview 

Galaxy clusters occupy a special position in the hierarchy of cosmic structure in many respects. 
Being the largest bound structures in the universe, they contain hundreds of galaxies and hot X- 
ray emitting gas and thus can be detected at large redshifts. Therefore, they appear to be ideal 
tools for studying large-scale structure, testing theories of structure formation and extracting 
invaluable cosmological information (cf. BohringerEJ, Schindlert3, Borgani & Guzzol). 

Below I will review a few issues related to cluster dynamics and the cluster large-scale 
environment. Using the APM cluster catalogue (Dalton et aZllHI) and methods calibrated using 
both optical and X-ray data I will present significant evidence that the cluster internal structure 
and dynamical state is strongly related to their large-scale distribution. 

This contribution is based on different works which are in preparation with collaborators 
of different parts of the project being S.Basilakos, S.Maurocordato, E.Slezak, C.Benoist and 
others. 



1.1 Cluster Internal Dynamics & Cosmology 

One of the interesting properties of galaxy clusters is the relation between their dynamical state 
and the underlying cosmology. Although the physics of cluster formation is complicated (cf. 
barazm , it is expected that in an open or a fiat with vacuum-energy contribution universe, 
clustering effectively freezes at high redshifts [z — Sl^ 1 — 1) and thus clusters today should 
appear more relaxed with weak or no indications of substructure. Instead, in a critical density 
model, such systems continue to form even today and should anpear to be dynamically active 
(cf. Richstone, Loeb & TurnerEl, Evrard et al. 0, Lacey & Cole0). Using the above theoretical 
expectations as a cosmological tool is hampered by two facts: 



Ambiguity in identifying cluster substructure: One has to deal with the issue of unambigu- 
ously identifying cluster substructure, since projection effects in the optical can conspire 
to make cluster images appear having multiple peaks/substructure. Alot of work has been 
devoted in attempts to find criteria and methods to identify cluster substructure (see ref- 
erences in Kolokotronis et aZEj). It is evident from all the available studies that there is 
neither agreement on the methods utilised nor on the exact frequency of clusters having 
substructure. 

Unknown physics of cluster merging: The clear-cut theoretical expectations regarding the 
fraction of clusters expected to be relaxed in different cosmological backgrounds break- 
down due to the complicated physics of cluster merging (cf. SarazinU) and especially due 
to the uncertainty of the post-merging relaxation time. In other words, identifying a cluster 
with significant substructure does not necessarily mean that this cluster is dynamically 
young, but could reflect an ancient merging that has not relaxed yet to an equilibrium 
configuration. 

Mohr et al. EH, Rizza et al. §1 and Kolokotronis et al. have investigated the frequency of 
cluster substructure using in a complementary fashion optical and X-ray data. The advantage 
of using X-ray data is that the X-ray emission is proportional to the square of the gas density 
(rather than just density in the optical) and thereforeit is centrally concentrated, a fact which 
minimises projection effects (cf. Sarazint3, Schindler a). The advantage of using optical data 
is the shear size of the available cluster catalogues and thus the statistical significance of the 
emanating results. Kolokotronis et aliB calibrated various substructure measures using APM 
and ROSAT data of 22 Abell clusters and found that in most cases using X-ray or optical 
data one can identify substructure unambiguously. Only in ~ 20% of the clusters they studied 
did they find projection effects in the optical that altered the X-ray definition of substructure. 
Their conclusion was that solely optical cluster imaging data can be used in order to identify 
the clusters that have significant substructure. 

However, as we discussed previously, it seems that in order to take advantage of the different 
rates of cluster evolution in the different cosmological backgrounds one needs (a) to find criteria 
of recent cluster merging and (b) calibrate the results using high-resolution cosmological hydro 
simulation, which will provide the expectations of the different cosmological models. 

Such criteria have been born out of numerical simulations (cf. Roettiger et al. and are 
based on the use of multiwavlength data, especially optical and X-ray data but radio as well (cf. 
Zabuldoff & ZaritskyB, Schindler 0, Sarazin, in this volume). The criteria are based on the 
fact that gas is collisional while galaxies are not and therefore during the merger of two clumps, 
containing galaxies and gas, we expect: (1) a difference in the spatial positions of the highest 
peak in the galaxy and gas distribution, (2) the X-ray emmiting gas, due to compression along 
the merging direction, to be elongated perpendicularly along this direction and (5) temperature 
gradients to develop due to the compression and subsequent shock heating of the gas. 

Figure 1 presents the smoothed optical (APM) and X-ray (ROSAT pointed observations) 
density distributions of 4 Abell clusters, out of which 3 show such evidence of recent merging, 
while A2580 seems to be in a relaxed state (for details of the smoothing procedure see Kolokotro- 
nis et a/E3). Note that BeppoSAX observations of Abell 3266 (de Grandi & MolendiH) has 
shown the existence of a temperature gradient dropping from 10 keV in the cluster core to 5 
keV at about 1.5 Mpc distance, consistent with a merging event, a fact which is also apparent, 
in figure 1, from the comparison of the optical and X-ray image of the cluster. 

1.2 Cluster Alignments and Formation Processes 

Another interesting observable, that was thought initially to provide strong constraints on the- 
ories of galaxy formation, is the tendency of clusters to be aligned with their nearest neighbour 





Figure 1: Optical APM (colour) and ROSAT Xray (contour) images of 4 Abell clusters. Peaks of the APM galaxy 
distribution is shown in blue. Note that appart from A2804, the rest are HRI images. 



as well as with other clusters that reside in the same supercluster (cf. Binggeliu, WestCil, Plionis 
@). Analytical work of Bondii in which clusters were identified as peaks of an initial Gaus- 
sian random field, has shown that such alignments, expected naturally to occur in "top-down" 
scenarios, are also found in hierarchical clustering models of structure formation like the CDM. 
These results were corroborated with the useof high-resolution N-body simulations by West et 
a/cll, Splinter et a0 and Onuora & Thomas c3. This fact has been explained as the result of an 
interesting property of Gaussian random fields that occurs for a wide range of initial conditions 
and which is the "cross-talk" between density fluctuations on different scales. This property is 
apparently also the cause of the observed filamentariness observed not only in " pancake" models 
but also in hierarchical models of structure formation; the strength of the effect, however, differs 
from model to model. 

There is strong evidence that the brightest galaxy (BCGs) in clusters are aligned with the 
orientation of their parent cluster and even with the orientationof the large-scale filamentary 
structure within which they are embedded (cf. StrubleE3, Westc3, Fuller, West & Bridges E-a). 
Furthermore, there is conflicting evidence regarding the alignment of cluster galaxies in general 
with the orientation of their parent cluster (cf. Djorgovskia, van Kampen & Rheelis, Trevese, 
Cirimele & FlinU). It may be that general galaxy alignments may be present in forming, 
dynamically young, clusters, while in relaxed clusters violent and other relaxation processes 




Figure 2: The distribution of galaxy 
position angles, within 0.75 h~ x 
Mpc, of A521. Position angles of 
galaxies with m < m* are shown in 



1 0.5 o -0.5 -1 red. The alignment with the clus- 

X ter position angle of 6 ~ 140° is ev- 

ident. 

may erase such alignment features. Such seems to be the case of the Abell 85/87/89 complex 
(see Durret et ali-H) and of Abell 521, a cluster at z ~ 0.25 which is forming at the intersection 
of two filaments (Arnaud et aln). Using wide-field CFHT imaging data of A521, Plionis et 
al {in preparation) have found statistical significant alignments not only of the predominantly 
bright but also of fainter galaxies with the major axis direction of the cluster (figure 2). It is 
interesting that this direction coincides with the direction of the nearest Abell cluster (see figure 
. Within the framework of hierarchical clustering, the anisotropic merger scenario of West 
, in which clusters form by accreting material along the filamentary structure within which 
they are embedded, provides an interesting explanation of such alignments as well as of the 
observed strong alignment of BCGs with their parent cluster orientation. Evidence supporting 
this scenario was presented in West, Jones & FormanE! in which they found, using Einstein data, 
that cluster substructures seem to be aligned with the orientation of their parent cluster and 
with the nearest-neighbouring cluster (see also Novikov et a/E3). Such effect has been observed 
also in numerical simulations of cluster formation for a variety of power-spectra (van Haarlem 
& van de WeygaertllH) 

Tidal effects between protoclusters could be an alternative explanation of galaxy and cluster 
alignments in the popular hierarchical clustering scenario. However, a back of the envelope 
calculation shows that in order to have a fractional effect / on a test particle at a distance of 
1 h^ 1 Mpc from a cluster of mass M c i, caused by a perturber at a distance R, one needs a 
perturber mass of: 

M p ~ fR 3 M c (1) 

Therefore for the typical intercluster nearest-neighbour separation of ~ 18 h^ 1 Mpc, one would 
get a 5% effect for M p ~ 300M C / ! This however does not exclude the possibility that tides, 
produced by the cluster itself, could effect the shape and orientation of cluster members (see 
Barnes Sz Efstathioui and Salvador-Sole & Solanesll for thorough studies of tidal effects). 
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Figure 3: The large scale environ- 
ment surrounding Abell 521. The 
major axis direction of A521 is 
pointing towards its nearest neigh- 
bour A517. 



Here I briefly present the methods used to determine the dynamical state of clusters, their shape, 
orientation and alignment. 



2.1 Substructure Measure - Cluster Centroid Shift 

Evrard et al. El and Mohr et al. @ have suggested as an indicator of cluster substructure the shift 
of the center-of-mass position as a function of density threshold above which it is estimated. 
The centroid- shift (sc) is defined as the distance between the cluster center-of-mass, (x ,y ), 
where x Q = J2 x iPi/J2 Ph Ho = UiPi/J2 A an d the highest cluster density-peak (x p ,y p ), ie., 



sc 



(x - x p ) 2 + (x Q - x p y 



(2) 



Notice here, that while the cluster center-of-mass changes as a function of density threshold, 
above which we define the cluster shape parameters, the position (x p ,y p ) remains unchanged. 
According to Kolokotronis et al 113, a large and significant value of sc is a clear indication of 
substructure. 

The significance of such centroid variations to the presence of background contamination 
and random density fluctuations are quantified using Monte Carlo cluster simulations in which, 
by construction, there is no substructure. For each cluster a series of simulated clusters is 
produced having the same number of observed galaxies as well as a random distribution of 
background galaxies, determined by the distance of the cluster and the APM selection function. 
The simulated galaxy distribution follows a King-like profile: 



S(r) oc 



1 + 



(3) 



where r c is the core radius. We use th 
and a determinations (cf. Girardi et al. 



ghted, by the sample size, mean of most recent r c 
), i.e., r c ~ 0.085 h^ 1 Mpc and a ~ 0.7. We do test 



the robustness of our results for a plausible range of these parameters (details can be found in 
Kolokotronis et a/lll). Naturally, we expect the simulated clusters to generate small sc's and in 
any case insignificant shifts. Therefore, for each optical cluster, 1000 such Monte-Carlo clusters 
are generated and we derive (sc) 

sim 

function of the same density thresholds as in the real 
cluster case. Then, within a search radius of 0.75 h~ l Mpc from the simulated highest cluster 
peak, we calculate the quantity: 

^ _ ( sc )o ~ ( sc )sim 

which is a measure of the significance of real centroid shifts as compared to the simulated, 
substructure-free clusters. Note that (sc) is the average, over three density thresholds, centroid 
shift for the real cluster. 



2.2 Cluster Shape Parameters 

To estimate the cluster parameters we use the familiar moments of inertia method (cf. Plionis, 
Barrow & FrenkB); with I n = £ Wi(rf - xf), I 22 = E Wi(rf - yf), I 12 = hi = -£ WiX^, 
where X{ and yi are the Cartesian coordinates of the galaxies and Wi is their weight. We, then 
diagonalize the inertia tensor solving the basic equation: 

det(hj - X 2 M 2 ) = 0, (5) 

where M 2 is the 2x2 unit matrix. The cluster ellipticity is given by e = 1 — where \ are 
the positive eigenvalues with (Ai > A2). 

This method can be applied to the data using either the discrete or smoothed distribution 
of galaxies (for details see Basilakos, Plionis & Maddoxi): 

• Method 1 (wi = 1): Galaxies within an initial small radius around the cluster center of 
mass are used to define the initial value of the cluster shape parameters. Then, the next 
nearest galaxy is added to the initial group and the shape is recalculated. Finally we 
obtain the cluster shape parameters as a function of cluster-centric distance. 

• Method 2 (wi = Sp/p): This method is based on applying the moment of inertia procedure 
on the smoothed galaxy density distribution which we obtain using a Gaussian or other 
kernel on a grid. 

Both methods have been found to provide consistent values of the cluster orientation but not so 
of the cluster ellipticity. 

In order to test whether there is any significant bias and tendency of the position angles to 
cluster around particular values we estimate their Fourier transform: 

/ 2 \ V2 

Cn= {N) J2 cos2n0 i (6) 

f 2 \ 1/2 

Sn = \ AT i z2 Sm 2nGi 
^ ' i=l 

If the cluster position angles, 0j, are uniformly distributed between 0° and 180°, then both C n 
and S n have zero mean and unit standard deviation. Therefore large values (> 2.5) indicate 
significant deviation from isotropy. 
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Figure 4: The fraction of clusters having substructure above the indicated significance level. 
2.3 Alignment Measures 

In order to investigate the alignment between cluster orientations, we define the relative position 
angle between cluster pairs by, (f>ij = |0j — OA. In an isotropic distribution we will have ((f>ij) — 
45°. A significant deviation from this would be an indication of an anisotropic distribution which 
can be quantified by (Struble &: Peebles H): 

N , 

^ = E^- 4 5 (8) 

i=l 

In an isotropic distribution we have (5) ~ 0, while the standard deviation is given by a = 
90/Vl2iV. A significantly negative value of S would indicate alignment and a positive misalign- 
ment. 

3 Results 

3.1 Cluster substructure 

Applying the above methodology to the ~ 900 APM clusters we find that about 30% of clusters 
have significant (> 3a) substructure. In figure 4 we present the differential and cumulative 
distribution of clusters as a function of substructure significance. Note that defining as having 
significant substructure those clusters with a > 2 or 2.5 increases the fraction to ~ 50% and 
40% respectively. Furthermore, changing the parameters of the Monte-Carlo clusters changes 
the actual cr-values, although their relative significance rank-order remains unaltered. Our 
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Figure 5: Upper panel: The distribution of relative position angles between nearest-neighbours. Lower panel: 

The distribution of APM cluster position angles. 



results are in accordance with a similar study of the BCS and REFLEX clusters in which a 
considerable number of clusters do show indications of significant substructures (Schuecker et 
al, in preparation). 



3.2 Cluster Alignments 

We have tested whether the well known nearest-neighbour alignment effect, present in the Abell 
clusters (cf. Bingelli 1982; Plionis 1994), is evident also in the poorer APM clusters. 

A necessary prerequisite for this type of analysis is that there is no orientation bias in the 
distribution of estimated cluster position angles. In the lower panel of figure 5 we present the 
corresponding distribution for the APM clusters. It is evident that the distribution is isotropic, 
as it is also quantified by the Fourier analysis. In the upper-panel of figure 5 we present the 
distribution of relative position angles, 5cj), between APM nearest-neighbours for two separation 
limits (one for all separations and one for separations < 10 hr 1 Mpc). It is evident that there 
is significant indication of cluster alignments in the small separation limit. 



0.35 



0.3 



i i r 



0.25 








A 
V 



10 



20 



30 




D ^20 h- 1 

cc 







4 

a 



6 



8 



Figure 6: Lower panel: Alignment signal as a function of substructure significance level. Blue dots correspond to 
nearest-neighbours and white dots to all neighbours within 20 h -1 Mpc. Upper panel: The mean cluster ellipticity 

as a function of a. 



3.3 Cluster Alignments versus Substructure 

We have correlated the alignment signal with the substructure significance indication in order 
to see whether there is any relation between the large-scale environment, in which the cluster 
distribution is embedded, and the internal cluster dynamics. 

In the lower panel of figure 6 we present the alignment signal, (<5), between cluster nearest- 
neighbours (blue dots) and between all pairs (open dots) with pair separations < 20 hr x Mpc. 
There is a strong correlation between the strength of the alignment signal and the substructure 
significance level. This result, based on the largest cluster sample available, supports the hierar- 
chical clustering scenario and in particular the formation of cluster by an isotropic merging along 
the filamentary structure within which they are embedded (cf. Westci El) . This is supported also 
by the fact that the mean ellipticity of clusters increases as a function substructure significance 
level, which is in agreement with the flatness being a result of accreting substructures along a 
preferred direction. 

If this view is correct then one would expect that clusters with significant substructure 
should be residing preferentially in high-density environments (superclusters), and this would 
then have an imprint in their spatial two-point correlation function. In the upper panel of figure 
7 we present the spatial 2-point correlation function of all APM clusters (open dots) and of those 
with substructure significance a > 4 (red dots). It is clear that the latter are significantly more 
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Figure 7: Upper panel: 2-point correlation function of all APM clusters (open dots) and of a > 4 clusters. Lower 
panel: The cluster correlation length and the mean APM cluster richness as a function of substructure significance. 

clustered. This can be seen also in the lower panel were we plot the correlation length, ro, as a 
function of a, which is clearly an increasing function of cluster substructure significance level. 
To test whether this effect could be due to the well-known richness dependence of the correlation 
strength, we also present the mean APM richness as a function of a and see that if any, there 
is only a small such richness trend. The conclusion of this correlation function analysis is that 
indeed the clusters showing evidence of dynamical activity reside in high-density environments, 
as anticipated from the alignment analysis. It is interesting that such environmental dependence 
has also been found for the cooling flow clusters with high mass accretion rates (Loken, Melott 
& Miller 0). 

3.4 Recent Evolution of Clusters ? 

An interesting exercise is to see whether any of the previously discussed features show an evo- 
lution with redshift, since in a low-density (flat or not) universe we do not expect to see recent 
significant evolution of the structural or dynamical parameters of clusters. As a first such test we 
present in figure 8 the redshift dependence of the ellipticity for those clusters that show evidence 
of significant substructure (a > 6). We find correlation coefficient R ~ 0.3 with probability 




of no correlation P < 10 -4 . A similar trend is found for the whole APM sample but of lower 
correlation and significance. It is interesting that Plionis et al. II found a similar trend for 
Abell clusters using the Lick galaxy catalogue. If this result survives a thorough investigation of 
possible systematic effects then it would be interesting to compare this observable with N-body 
simulations of different cosmological models in an attempt to constrain the models. 

4 Conclusions 

We have presented evidence, based on the large APM cluster sample, that there is a strong link 
between the dynamical state of clusters and their large-scale environment. Clusters showing 
evidence of dynamical youth are significantly more aligned with their nearest neighbours and 
they are also much more spatially clustered. This supports the hierarchical clustering models 
and in particular the anisotropic merger scenario of West c3. 
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